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Gradual localization of charge carriers was studied in a series of micro-size samples of monolayer 
graphene fabricated on the common large scale film and irradiated by different doses of C"*" ions with 
energy 35 keV. Measurements of the temperature dependence of conductivity and magnetoresistance 
in fields up to 4 T showed that at low disorder, the samples are in the regime of weak localization and 
antilocalization. Further increase of disorder leads to strong localization regime, when conductivity 
is described by the variable-range-hopping (VRH) mechanism. A crossover from the Mott regime to 
the Efros-Shklovskii regime of VRH is observed with decreasing temperature. Theoretical analysis 
of conductivity in both regimes showed a remarkably good agreement with experimental data. 

PACS numbers: 73.22.Pr 


INTRODUCTION 

Graphene, a sheet of sp^ -bonded carbon atoms, for 
several years already continues to be in the focus of atten¬ 
tion of physics community, due to its potentially transfor¬ 
mative impact across a wide range of applications includ¬ 
ing advanced electronics and sensing. Carrier scattering 
in graphene can be due to all kinds of disorder, including 
ripples in the graphene layer, point defects and their asso¬ 
ciated short-range potentials, charged impurities residing 
in the supporting substrate, and adsorbed atoms on the 
surface. Investigation of the influence of disorder on the 
properties of graphene is attracting a tremendous inter¬ 
est due to possibility to modify this novel and promising 
material using weak or strong localization of charge car¬ 
riers. By controllably introducing defects into graphene, 
one may be able to understand how these mechanisms 
limit transport. 

Previosly, there were observed separately either weak 
localization (WL) or different kinds of Variable-Range- 
Hopping (VRH) conductivity of strongly localized carri¬ 
ers in graphene samples disordered by different methods 
like doping, oxidation, ion irradiation (see, for example, 
[IH). However, we are not aware of observations of all 
regimes of localization with gradual increase of disorder¬ 
ing in graphene. In this paper we report the results of 
study of the localization process in monolayer graphene 
(MG) samples subjected by different doses of ion irradi¬ 
ation. 

The initial large size (5x5 min) specimens were sup¬ 
plied by ” Graphenea” company. Monolayer graphene was 
produced by GVD on copper catalyst and transferred to 
a 300 nm Si02/Si substrate using wet transfer process. It 
is specified in the certificate, that the sheet resistivity of 
the specimen is 350 Ohm/sq. On one specimen, gold elec¬ 


trical contacts were deposited directly on the graphene 
surface. This sample is marked as 0. The resistivity of 
the sample 0 was 380 Ohm/sq which is close to the data 
in the certificate. On the surface of the other specimen, 
six groups of mini-samples (200 x 200 /rm) were fabri¬ 
cated by means of electron-beam lithography (EBL) as 
well as electrical contacts (5 nm Ti and 45 nm Pd) for 2- 
probe measurements. The samples from the first group, 
marked as sample 1 were not irradiated, while 5 other 
groups were subjected to different doses (from 5 x 10^3 
up to 1X10^® cm“^) of irradiation by G’*' ions with energy 
35 keV. Ion irradiation was performed on HVEE-350 Im- 
planter, ions were obtained from the hollow cathode 
ion source by CO 2 decomposition. 

The choice of the irradiation conditions was not oc¬ 
casional. They are the same as in the Ref. Q where 
the Hall effect was measured in monolayer graphene sub¬ 
jected to ion irradiation. Resistances of the initial sample 
in Ref. Q and our initial sample 0 are the same. It was 
shown in that both Hall coefficient and mobility are 
temperature independent up to high dose of irradiation. 
Moreover, it was shown that the sheet Hall concentration 
of charge carriers (about 10^^ cm“^) does not depend at 
all on the dose of irradiation. We were based on these 
results in our decision to prefer the 2-probe geometry in¬ 
stead of the Hall bar geometry, because this allows us 
to measure more samples on the same sample holder. 
As a result, we were able to enhance the reproducibil¬ 
ity of the obtained results. Comparison of the graphene 
resistance measured W 2-probe method and by 4-probe 
method presented in [6| showed that in this system, con¬ 
tact resistance is insignificant. 

In our previous work Q, a concentration of structural 
defects Nd was determined for each group of samples us¬ 
ing measurements of the Raman scattering. These values 
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are shown in the inset in Fig. [TJ (It turns out that sample 
1 is also slightly disordered due to EBL process). Mea¬ 
surements of the current-voltage characteristics {I — V ) 
for all samples, performed at room temperature, showed 
that for samples 5 and 6, 1—V is strongly non-linear even 
at very small current. That is why in this paper, the tem¬ 
perature dependences of resistance R(T) are shown only 
for samples 0-4. 

The resistance was measured by two-probe method in 
helium cryostat down to 1.8 K in zero magnetic field and 
in magnetic fields up to 4 Tesla. Fig. [ijshows the general 
picture R{T) for all samples. 

The sample 0 shows typical metallic behavior, when 
R slightly decreases with decrease of T. For sample 1 
R slightly increases with decreasing T, which is charac¬ 
teristic for ’’dirty” metals. For other samples R changes 
with T exponentially, which is characteristic for strongly 
localized carriers. 
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FIG. 1: Resistivity of disordered monolayer graphene sam¬ 
ples as a function of temperature. Inset shows the density of 
structural defects in samples [^. 


WEAK LOCALIZATION 

Fig. shows the experimental dependences of mag¬ 
netoconductance (MC) of sample 1 in wide temperature 
interval, from 300 K down to 1.8 K. Plot of the temepa- 
ture dependence of conductivity on the scale a vs. InT 
(Fig. [3]) shows the logarithmic temperature behavior of 
cr at low T, characteristic for regime of WL 0 , with 
tendency to saturation at very low temperatures. 

WL regime of conductivity in monolayer graphene has 
important features due to the facts that charge carriers 
are chiral Dirac fermions, which are reside in two inequiv¬ 
alent valleys at the K and K' points of the Brillouin 
zone. Due to chirality, Dirac fermion acquires a phase 
of TT upon intravalley scattering, which leads to destruc¬ 
tive interference with its time-reversed counterpart and 
weak antilocalization (WAL). Intervalley scattering leads 
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FIG. 2: Magnetoconductance of sample 1 as function of mag¬ 
netic field; solid lines - experiment, dashed lines - formula © 
with htted parameters. 


to restoration of WL because fermions in K and K' val¬ 
leys have opposite chiralities. 

Quantum corrections to the conductivi ty of g raphene 
have been intensely studied theoretically |lll4l9l |. It was 
predicted that at relatively high temperatures WAL cor¬ 
rections will dominate, while with decreasing T the WL 
corrections will dominate. There were several experimen¬ 
tal papers reporting logarithmic dependence of conduc¬ 
tivity on tenmerature and magnetic field at low temper¬ 
atures However, in our sample I, the logarith¬ 

mic dependence is observed in wide temperature interval, 
starting from 300 K, which gives an opportunity to check 
in a very detailed way the theoretical predictions. 

For the MC the theory 0 predicts 
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where "0 is the digamma function, is the coherence 
time, T~^ is the intervalley scattering rate, is the 
combined scattering rate of intravalley and intervalley 
scattering and of trigonal warping. 

Fitting Eq. to experimental data for inagnetocon- 
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FIG. 3: Conductivity of sample 1 as function of tempera¬ 
ture. Circles present experimental data, squares present Eq. 
(H) with parameters determined from fitting the magnetocon¬ 
ductance {A was chosen to be 11.7 • c^/tt/i). The inset: The 
values of (solid line) and Bi (short-dashed line) and B* 
(long-dashed line, right axis). 


ductance of sample 1 at different temperatures is illus¬ 
trated on Fig. [21 In the process of fitting we were able 
to extract all three parameters entering the equa¬ 

tion. These parameters shown in inset in Fig. [31 It turns 
out, that they are temperature-dependent which was not 
predicted by theory. Using these parameters we can cal¬ 
culate conductance at zero magnetic field according to 
Eq. (10) of Ref. 13 , which can be rewritten in the form 
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where A is a constant, dependent upon the unit of mag¬ 
netic field (chosen as 1 T). We compare Eq. ([2|) with the 
experimentally measured conductivity a{B = 0,T). The 
comparison is presented on Fig. [31 The good agreement 
proves the correctness of the chosen parameters. 

The inset to Fig. [3] shows that, apart from the low¬ 


est temperatures, B^, 
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T. Mechanisms that 


can give the dependence t,„ ~ l/T are: electron-electron 
scattering in dirty limit |10j , electron-phonon scatte ring , 
[ 2 ^ 21 1, and electron-flexural phonon interaction (22l [. 


The saturation of Tip at low temperatures is well known in 
classical 2d systems and may be connected with existence 
of dephasing centers (for example, magnetic impurities) 

^m- 

Obtained value of Btp allow us to determine the values 
of dephasing length L,p = Drip = ^Jhf^B^. When the 
temperature decreases from 300 K to 3 K, the dephasing 


length Lip increases from 7 nm to 70 nm and then sat¬ 
urates. The maximal value of Lip allows us to estimate 
the density of dephasing centers as 2 x 10^° cm“^. 


STRONG LOCALIZATION 


Lets discuss now samples 2-4 with pronounced insulat¬ 
ing behaviour. Plotting the data on the Arrhenius scale 
Ini? vs. l/T showed that energy of activation continu¬ 
ously decreases with decreasing T which is characteristic 
for the variable-range-hopping (VRH) conductivity (25| . 
There are two kinds of VRH depending on the structure 
of the density-of-states (DOS) g{e) in the vicinity of the 
Fermi level (FL) /r: when g{e) = g{fi) =const, R{T) is 
described by the Mott (T~^/^) law in the case of two- 
dimensional (2d) conductivity: 

i?(r) =i?oexp(rM/T)i/3, TM = CM[gi^^i)aT". (3) 

Here Cm = 13.8 is the numerical coefficient [ 2 ^, a is the 
radius of localization. 

The Coulomb interaction between localized carriers 
leads to appearance of the soft Coulomb gap in the vicin¬ 
ity of FL which in the case of 2d has a linear form 

5(e)~|e-/7|(eV«)-^ (4) 


where k is the dielectric constant of the material. This 
leads to the Efros-Shklovskii (ES) VRH or law: 

i?(r) = i?oexp(r£;s/T)i/2, TEs = CEs{e^l^^a), (5) 


where the numerical coefficient Ces = 2.8 |25l |. 

Coulomb interaction can alter the DOS only near the 
FL. Far from FL, the DOS is restored to its initial value, 
which is approximately equal to g{g), see inset in Fig. (4] 
Denoting the half-width of the Coulomb gap as A one 
can conclude, therefore, that T <g; A, ES law has to be 
observed, while in the opposite case (T A), the Mott 
law should dominate. 

There is a number of reports about observation of ei¬ 
ther Mott or ES laws in different disordered graphene- 
based materials [H-iii. We show that in samples 3 and 
4, both VRH laws are observable at different tempera¬ 
tures. (For sample 2, the VRH regime will be observed 
at lower temperatures). In Figs. |4]and|5l logi? is plotted 
versus and At high temperatures, depen¬ 
dences i?(T) are straightened on the scale while 

at low temperatures they are straightened on the scale 
T- 1 / 2 ^ The latter shows the approach to the ES law 
which should be observed at the lowest temperatures. 
These plots allow us to determine both parameters Tm 
and Tes (Table 1) and calculate the temperature Tc of 
deviation from law to law in the case of 2d 

conductiviW similarly to the calculation of Tc for 3d con¬ 
ductivity [30|. In VRH, only localized states in an opti- 
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S # 

Tm(K) 

Tes{K) 

Tc(K) 

A(K) 

3 

308 ± 54 

50 ±7 

11.4 ±0.8 

12 ± 1.2 

4 

5962 ± 305 

490 ± 12 

29 ±5 

60 ±6 


TABLE I: Hopping conductivity parameters for samples 3,4. 

mal band of width e(T) near the Fermi level are involved 
in the hopping process. The band becomes continuously 
narrower with decreasing T. Hopping resistance is deter¬ 
mined by the critical parameter 

R = Roexp^c, = i^r/a) + {e/T). (6) 

Here energy and temperature are measured in the same 
units, r is the mean distance of hopping. In the Mott 
regime, g{e) = g{g) = const and, therefore, the total 
number of states in the optimal band is N{T) = g{g)e 
and the mean distance between states in two dimensions 
is r ss [g{g)t]~'^/‘^. Substituting in dH]), one can find e(T) 
from the minimal value of determined from d^c/de = 0: 


T(K) 



FIG. 4: logR for samples 2-4 plotted versus T . 


e(T) = r2/3[g(/,)a2]-i/3. (7) 

This gives the relationship between T and the width of 
the optimal band: T = [g{fj,)a^]^^'^€^^^. At the crossover 
temperature Tc, e = A, so 

T, = [5(/i)a2]i/2A3/2. (8) 


Being inside the Coulomb gap, the crossover temperature 
can be determined from g{e = A) = g{g.) which gives 
A = g[g){e^ / kY . Substituting into Eq. dS]) and using 
expressions for Tm and Tes, Eqs. © and ®, we get 

Tc = {ClJCls)(Tls/Tli) « mns/Tlt) (9) 


The functional proportionality of Tc to the ratio 
{T^s/'^m) been obtained earlier in Ref. 31, 3^, but 
with significantly different numerical coefficient. The lat¬ 
ter is, however, crucial for comparison with experiment. 
The values of Tc calculated from Eq. m for samples 3 
and 4 are given in Table 1 and shown as arrows in Fig. 
m The good agreement shows the correctness of the ob¬ 
tained numerical coefficient. We can also estimate the 
width of the Coulomb gap. In calculation, the value of 
K = 2.45 for the monolayer graphene on the Si02 surface 
was used as k = (ki -I- K 2 ) j ‘ 2 , where ki = 3.9 (for Si02) 
and K 2 = 1 (for air). One can see that indeed, the ES 
law is observed when T < A. 

Comparison of samples 3 and 4 shows that increase 
of the density of defects Ne leads to stronger localiza¬ 
tion which manifests in a significant decrease of a and 
increase of the energy band needed for hopping. This 
looks like increase of the amplitude of the random po¬ 
tential relief in classical semiconductors induced by ran¬ 
domly distributed positively charged donors and nega¬ 
tively charged acceptors in the case of compensation. 
Hence we assume that structural defects in graphene are 
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FIG. 5: logR for samples 2-4 plotted versus T 


of amphoteric impurity action [33j, i.e. they can be either 
acceptors or donors and compensate each other. Thus, 
increase of Ne leads to the random potential relief ampli¬ 
tude increase. One can also assume that this phenomena 
may be connected with the fact that some point defects 
(say, vacancies) may produce complex associative cen¬ 
ters with other defects. These associations can show an 
amphoteric impurity action, whereas the individual com¬ 
ponents are not amphoteric 34l |. 

In conclusion, a gradual transformation of conductiv¬ 
ity measured in a wide interval of temperatures (300 - 1.8 
K) and magnetic fields (up to 4 T) was observed in a se¬ 
ries of monolayer graphene samples subjected by different 
dose of ion irradiation. Increasing the density of struc¬ 
tural defects Ne induced by irradiation has led to change 
the mechanism of electron transport from metallic con¬ 
ductivity to the regime of WL and finally to the VRH of 
strongly localized carriers. It is shown that WL regime in 
slightly disordered sample starts from high (room) tern- 










































5 


peratures. 

Comparison of experimental magnetoconductance 
curves with theory allowed us to find the parameters 
which determine the logarithmic temperature depen¬ 
dence of conductivity in WL regime. In VRH, a crossover 
from the Mott law to the Efros-Shklovskii law was ob¬ 
served in the same samples with decreasing tempera¬ 
ture. The calculated crossover temperatures are in good 
agreement with experimental values. It is suggested that 
strengthening of localization with increase of No can 
be explained by amphoteric impurity action of graphene 
structural defects induced by ion irradiation. 
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